PIXELATIONS OF PLANAR SEMIALGEBRAIC SETS 



LIVIU I. NICOLAESCU AND BRANDON ROWEKAMP 

Abstract. We describe an algorithm thiat associates to each positive real number r and each finite 
collection Cr of planar pixels of size r a planar piecewise linear set Sr with the following additional 
property: if Cr is the collection of pixels of size r that touch a given compact semialgebraic set 5", then 
the normal cycle of Sr converges to the normal cycle of 5* in the sense of currents. In particular, in the 
limit we can recover the homotopy type of 5* and its geometric invariants such as area, perimeter and 
curvature measures. At its core, this algorithm is a discretization of stratified Morse theory. 
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Introduction 

This paper is a natural sequel of the investigation begun by the second author in his dissertation 
[12]. To formulate the main problem discussed in [12] and in this paper we need to introduce a bit of 
terminology. 

For e > we define an e-pixel to be a square of the form 

[(m - me] x [(n - ne] C M^, 

where m,n £ Z. The number e is called the resolution. A pixelation is a union of finitely many pixels. 
The e-pixelation of a set 5 C is the union of all the pixels that touch S. We denote it by Pe{S). 

The Main Problem. Produce an algorithm that associates to Pe{S) a PL-set Ss which approximates 
S very well as e \ 0. More precisely, for e > sufficiently small, the approximation Sg must have 
the same homotopy type as S and the curvature features of Sg must closely resemble those of S. 

The pixelation Pe{S) is a PL-set. It is however very jagged and there is no hope that its curvature 
properties are similar to those of S. In fact there is a more insidious reason why the pixelation is a poor 
approximation for S. 

Consider the pixelation of an angle with vertex at the origin whose edges have slopes | and 1. Figure 
1 shows that this pixelation is not contractible and in fact its first Betti number is 2. These two "holes" 
won't disappear at any resolution because all the pixelations of this angle are rescalings of each other. 
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Figure 1 . The pixelation of the angle at the origin formed by a line of slope 1 and a 
line of slope |. 

In [12] the second author solved the Main Problem in the special case when S itself is a PL-set. The 
resulting algorithm is based on two key principles. 

Principle 1. If S is the graph of a piecewise C^-function /, then the pixelation of the graph of / has 
the same homotopy type as the graph, i.e., it is contractible. Moreover, every column of the pixelation 
is connected, where by column we understand the intersection of the pixelation with a vertical strip of 
the form {{m — l)e < x < me}. We fix a function a : (0, oo) — )• Z>o, called the spread, such that 

lim ea(e) = and lim ecj(e)^ = oo (cr) 

For every e > we obtain by linear interpolation a PL function (resp. /J^) whose graph is produced 
by connecting with straight line segments the centers of the top (resp. bottom) pixels of every a{e)-th 
column of the pixelation of the graph of /. 



Figure 2. Linear interpolations with spread a = 3. 

The result of the algorithm is the PL-region between the graphs of /?£ and r^. This is a very 
narrow two dimensional PL set very close to the graph of /. Moreover, the condition (a) gurantees 
that the curvature of resembles that of 5. 

An identical strategy works when 5 is a set of the form 

S={{x,y)eM.^; xe[a,b], ; P{x) < y < t{x)} , 

where /3,r : [a,b] — M are Lipschitz continuous, piecewise C^-functions such that /3(x) < t(x), 
Vx G [a, b]. 
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We will refer to these two types of sets as elementary. Thus, the Main Problem has a solution for 
elementary sets. 

Principle 2. Consider the linear map h{x, y) = x. The Morse theoretic properties of the restriction of 
h to Pe{S) closely mimic the Morse theoretic properties of the restriction of /i to 5 if e is sufficiently 
small. Here are the details. 

For G M denote by ns{xo) the number of connected components of the intersection of S with 
the vertical line {x = xq} and denote by 3s the set of discontinuities of the function x i— )• ns{x). Then 
3s is a finite subset of M and there exists 7 > such that for any r G (0, 7) the set S'^. obtained from S 
by removing the vertical strips {\x — j\ < r}, j G 3s^ is a disjoint union of elementary regions. 

The set ^Js is difficult to determine from a pixelation, but one can algorithmically produce a very 
small region containing it. Here is roughly the strategy. 

For e > and xq G M we denote by ns^e{xo) the number of connected components of the intersec- 
tion of the vertical line {x = xq} with the pixelation Pe{S). We denote by 3s,£ the set of discontinuities 
of the function x 1— )• ns e(x). The set 3s,£ is finite and one can prove the following fact. 
(D) There exist h > 0, uq > and kq G (0, 1], depending only on S such that for e < ^ we have 

dist{3 s,ds,e) < iyoe^°. 
We define the noise region to be the set 

:= {x G M; dist{x,3s,e) < 2i^e''° }. 
For £ sufficiently small, the noise region is a finite union of disjoint compact intervals 

Ij{e), j = l,...,N ■.= #3s, 

called noise intervals. 

We denote by P^{S) the closure of the set obtained from Pe{S) by removing the vertical strips 
{x G /j(e)}, j = 1,. . . ,N. Each of the connected components of P^{S) is the pixelation of an 
elementary set and as such it can be PL-approximated using Principle 1. 

The approximation of the noisy regions, i.e., the intersection of Pe{S) with the above vertical strips 
is rather coarse. Every component of such a region is approximated by the smallest rectangle that 
contains it. Here by rectangle we mean a region of the form [a, b] x [c, d], a < b, c < d. 

It turns out that the approximation Se of S obtained in this fashion from Pe{S) is very good in the 
sense that the normal cycle of Se converges in the sense of cun^ents to the normal cycle of 5. For a nice 
introduction to the subject of normal cycles we refer to [9]. A brief description of this concept can also 
be found on page 20 of this paper. 

The goal of this paper is to extend the above program to the more general case of compact, semi- 
algebraic subsets of M?. While Principle 1 extends with little difficulty to the semi-algebraic case, 
Principle 2 requires a more delicate analysis. This requires that 5" be a generic semialgebraic set in the 
sense that the restriction to S of the linear function h{x, y) = xbe a. stratified Morse function, [4, 1 1]. 
Principle 2 is proved using results in real algebraic geometry. 

The PL approximation S'e of S is obtained as before, using the two principles. To prove that the nor- 
mal cycle of Se converges in the sense of currents to the normal cycle of S we rely on an approximation 
theorem of J. Fu, [5]. That theorem states that the convergence of the normal cycles is guaranteed once 
we prove two things. 

• Uniform bounds for the perimeter and total curvature of Se- 

• For almost any closed half -plane H we have 

iimx{HnSe)=x{Hns), 

where x denotes the Euler characteristic. 
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Of the above two facts, the second is by far the most deUcate, and its proof takes up the bulk of this 
paper. 

Let us say a few words about the organization of the paper. In Section 1 we introduce the terminology 
used throughout the paper. Principle 1 is proved in Section 2, while Principle 2 is proved in Section 
3. Section 4 is devoted to the proof of the convergence of normal cycle. The paper also contains two 
appendices. In Appendix A we collect a few basic facts of real algebraic geometry usd in this paper 
together with a few other technical results. In Appendix B we give a more formal description of the 
approximation algorithm in a way that makes it easily implementable on a computer. The present paper 
can be read independently of [ 1 2] which served both as motivation and as inspiration. 

1. Basic facts 
We begin by recalling some basic notions introduced in [12]. 
Definition 1.1. (a) Let e > 0. Then we define an e-pixel to be the square in of the form 

Sijie) = [(i - l)e,ie)] x [(j - l)eje)] C R\ ij G Z. 

(b) A union of finitely many e-pixels is called an e-pixelation. The variable e is called the resolution of 
the pixelation. 

(c) For any compact subset 5" C we define the e-pixelation of S to be the union of all the e-pixels 
that intersect S. We denote the e-pixelation of S by Pe{S). The pixelation of a function / is defined to 
be the pixelation of its graph r(/). We will denote this pixelation by Pe(/). □ 

Observe that if || • ||oo : ^ [0, oo) is the norm 

\\{x,y)\\ := max{|x|, |y|}, 
then the pixel Sij{e) can be identified with the || • ||oo-closed ball of radius e/2 and center 

Therefore a pixelation of the compact set S can be thought of as set of points chosen neai^ S. 

Definition 1.2. Fix e > and a compact set S* C M^. 

(1) A point a G M will be called e-generic if x G M \ eZ. For such a point a we denote by /^(a) 
the interval of the form ( ne, {n + n G Z that contains a. 

(2) For a < bwe define the vertical strip 

§a,b ■■= ia,b) X R 
For every A: G Z we denote by ^ the vertical strip 

( te, (/c + l)e ) XR = §'ke,{k+l)e 

For any e-generic point a G M we denote by Se(a) the strip 

/,(a) X R = k:= [x/k\. 

(3) A column of Pe{S) is the intersection of Ps{S) with a vertical strip fc, k £ Z. The connected 
components of a column are called stacks. 

(4) For every e-generic a G R, we define the column of a pixelation Pe{S) over a to be the set 

CsiS,a) ■.= §,ia)nPsiS). 

In other words, Ci;{S, a) is the union of the pixels in PsiS) which intersect the vertical line 
{x = a}. When S is the graph of a function /, we will use the notation Ce(/, a) to denote the 
column over a of the pixelation Psi f )- 
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(5) We define the top, bottom and respectively height of a column to be the quantities 

Te{S,a) := sup{y : (a,y) G Ce{S,a)}, 
Be{S,a) := mf{y : {a,y) G CeiS,a)}, 

and respectively 

K{S,a) := ^{Te{S,a)-Be{S,a)). 

As in the other definitions when S is the graph of a function / we will use the notations 

Te{f,a),B,{S,a) and he{S, a). 

□ 

We have the following result,[12, Thm. 2.2] 

Theorem 1.3. If f : [a,b] ^ M. is a continuous function, then for any e > the columns of the 
e-pixelation of the graph of f consist of single stacks. □ 

In this paper we will be concerned with pixelations of generic planar semialgebraic sets, where the 
genericity has a very precise meaning. To describe it we need to introduce some terminology from 
stratified Morse theory, [4, 1 I]. 

For any subset X C we denote by cl{X) its closure and by dtopX its topological boundary, 

atop(x) ■.= ci{x)\x. 

We define a good stratification of a compact semialgebraic set 5" C to be an increasing filtration 

3": C C = 5 

satisfying the following properties. 

• Each of the sets F^'^\ i = 0, 1, 2 is closed. 

• dimF^*) < i = 0, 1, 2. In particular F^^'^ is a finite collection of points called vertices of the 
good stratification. 

• The connected components of F^^^ \ F^^^ ai^e open real analytic arcs, i.e., images of injective 
real analytic maps (0, 1) — > M^. We will refer to these components as the arcs or the edges of 
the stratification. 

• The connected components of F^^^ \ F^^' are open subsets of R^. They are called the/ace^ of 
the stratification. 

atop ( F^^) \F^'))<Z fW , Stop \ ) C F W . 
If t; is a vertex of a good stratification of a compact semialgebraic set we define Coo(i', S) to be the 
set of one dimensional subspaces of C such that there exists an arc A of the stratification with 
the following properties. 

• V e cl{A). 

• There exists a sequence of points Vn & A such that as ?i — )• oo we have Vn ^ v and the tangent 
spaces Tv„A converge to Lqo- 

Suppose that S C is a compact semialgebraic set equipped with a good stratification 3", and 
/ : — )• M. A point p G 5 is said to be a critical point of the restriction /|s if either p is a vertex, 
or p is the critical point of the restriction of / to an arc or to a face. The critical point p is said to be 
nondegerate if it satisfies one of the following conditions. 
(Co) The point p is a vertex and for any Lqo G Coo{p, S), the differential of / at p does not vanish 
along Loo- 
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(Ci) The point p belongs to an arc A of the stratification and as such it is a nondegenerate point of 

fU 

(C2) The point p belongs to a face F of the stratification and as such it is a nondegenerate point of 
/If. 

A function / : — > M is said to be a a stratified Morse function with respect to the semialgebaric 
set S equipped with the good stratification 9" if all its critical points are nondegenerate, and no two 
critical points lie on the same level set of /. 

A compact semialgebraic set 5 C is called generic if it admits a good stratification 9" such that 
projection onto the x-axis {x, y) 1— )• 2; is a stratified Morse function with respect to {S, 9"). Denote this 
projection by h. 

Observe that if 3" is a good stratification of a compact semialgebraic set 5 C M^, then p is a critical 
point of h relative to {S, 3") if either p is a vertex of the stratification, or p is a point on an arc of 3" 
where the the tangent space is vertical. 



A 




Figure 3. The curve on the right is generic, while the curve on the left is not. 

In Figure 3 we have depicted two one-dimensional planar semi algebraic curves (arcs of circles). 
The marked points are critical points of h. The point A on the left-hand side curve is a degenerate 
critical point because the condition (Co) is violated: the vertical line is tangent to the curve at that 
point. 

In the left-hand side of Figure 4 we have depicted further examples of pathologies prohibited by the 
genericity condition. (The pathologies involve the points with vertical tangencies.) The right-hand side 
depicts generic sets obtained by small perturbations from the nongeneric sets in the left-hand side. 

A A 



Figure 4. The curves on the left-hand side are nongeneric. They become generic 
after a small perturbation. 



2. Approximations of elementary sets 

In this section we study the pixelations of simple two dimensional sets. 

Definition 2.1. Fix e > and a compact set S. 

(1) An e-profile of 5 is a set Ilg of points in the plane with the following properties. 

(a) Each point in Ilg is the center of an e-pixel that intersects S. 

(b) Every column of Pe{S) contains precisely one point of lie. 
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(2) The top/bottom e-profile is the profile consisting of the centers of the highest/lowest pixels in 
each column of Pe{S). 

(3) An e-sample of 5 is a subset of an e-profile. An upper/lower e-sample of S is an e-sample of 
the upper/lower e-profile of S. 

□ 



Definition 2.2. Suppose pi, . . . ,pn is a finite sequence of points in M^. (The points need not be 
distinct). We denote by 

{pi,P2,. ■ ■ ,Pn) 

the PL curve defined as the union of the straight line segments [pi,p2]i • • • i [Pn-i,Pn]- □ 

Observe that each e-profile He of a set is equipped with a linear order <. More precisely if pi,p2 
are points in He, then 

Pi ■< P2 <^=^ X{pi) < X{p2), 

where x : — ^ M denotes the projection (x, y) ^ x. In particular, this shows that any e-sample of S 
carries a natural total order. 

Definition 2.3. If H is an e-sample of S, then the PL-interpolation determined by the sample H is the 
continuous, piecewise linear function L = obtained as follows. 

• Arrange the points in E in increasing order, with respect to the above total order, 

^ = {^0 ^6 ^6 -< ••• ^Cn}, n + l = \E\. 

• The graph of L= is the PL-curve {£,o,£,i, ■ ■ ■ , Cn)- 

□ 



In applications, the sample sets E will be chosen to satisfy certain regularity. 

Definition 2.4. (1) A spread function is an increasing function (0,oo) 9 e i— )• a{e) G Z>o with 
the following properties. 

lim a{e) = oo, (2.1a) 
limea(e) = 0. (2.1b) 
(2) If o" is a positive integer and Ilg is an e-profile, then an e-sample with spread a is a subset 

H = Uo-<----<Ucn,(s) 

such that the following hold. 

• The points and ^„ are the left and rightmost points in the profile. (That is for each 

p £ He, x{^o) < X{p) < X{in)-) 

• For any p G He, there exists ^ G H such that \x{p) — x{^)\ < ea. 
• 

-ea < \x{Ck) - x{Ck-i)\ < eo-, V/c = 1, . . . , n. 

□ 



Definition 2.5. A subset 5 C is said to be elementary (with respect to the x-axis) if it can be defined 
as 

5 = 5(/3,r) := {(x,2/) :xG [a,6],/3(x) <x<r(x)}, 
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where /3,r : [a,b] — >■ M ai^e continuous semialgebraic functions such that < t{x), Vx G [a,b]. 
The function /3 is called the bottom of S while r is called the top of S. If 

/3(x) < t{x), Vx G (a, 6), (2.2) 

then the elementary set is said to be nondegenerate. If 

/3(x) =r(x), Vx G [a, 6], (2.3) 

then the set S is called degenerate. The elementary set is called mixed if both sets 

{xG(a,6); r(x) - /3(x) > } and {xG(a,6); t(x) - /3(x) = } 

are nonempty. □ 

Observe that an elementary set S'(/3, r), /3, r : [a, 6] — M admits good partitions, i.e., partitions 

a = Co < ci < ci • • • < Cn = 6, n > 2 

such that each of the elementary sets [cj_i, q] x M n S{P, r) is either degenerate or nondegenerate. 
The good partition with minimal n is called the minimal good partition; see Figure 5. 





T 


1 
1 

1 


















> > 



Figure 5 . The minimal good partition of a mixed elementary set. 

In the remainder of this section S will indicate an elementary set. We first note that like the pixelation 
of a function, each column of P£{S) contains only one stack. 

Proposition 2.6. If S = S'(/3, r) is an elementary set, then for every x G [a, b] \ eL, the column 
Cs{S,x) consists of exactly one stack. 

Proof. Fix an e-generic x G [a,b]. By Theorem 1.3 the columns Ce{P,x) and Cs{t,x) consist of 
single stacks. If these two columns intersect, then the conclusion is obvious. If they do not intersect, 
then any pixel in the strip Se{x) situated below the stack Ce(r, x) and above the stack Cs{f3, x) is a 
pixel of P£{S). This proves that the column Ce{S, x) consists of a single stack. □ 

Definition 2.7. Fix e > and an elementary set S = S{f3,T). Suppose that Ef are compatible 
upper/lower samples of S 

The PL-approximation of S determined by these two samples is the compact PL-set bounded by the 
closed PL-curve 



□ 
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The total curvature of a immersion 7 : [a, 6] — )• is defined as follows. Define a -function 
6* : [a, 5] M such that 

^7W = (cos(0(t),sin(0(t))), 
where a dot " denotes the t-derivative. We set 

K{t) := \m. 

The scalar K{t) is called the curvature of 7 at the point j{t). We define the total curvature of 7 to be 

K(7) := / (it. 

J a 

Suppose now that 7 : [a, b] — )■ continuous and piecewise C^-immersion, i.e., there exist a finite 
subset {to, . . . , tj/} C [a, 6], such that 

a = to < *i < • • • < tu~i <ty = h, 

and the restriction 7j := '^\[[t^_-^^t^\ is a C^-immersion for any i = 1, . . . ,v. The curvature of 7 at a 
jump point 7(ti) is the quantity 

<ti) = mt) - 9{tT)\ := I liin e{t) - lira e{t)\. 
We define the total curvature of 7 to be 

1=1 i=l 

'0, 7(b) # 7(a) 



+ 



l^(io+)-^(i^)l, 7(b) = 7(a). 



For more details we refer to [8] and [9, §2.2]. 

We define a semialgebraic arc to be the image of a continuous, injective semialgebraic map 

V? : [a, 6] 

Suppose that 93 : [a, b] — is a continuous, injective, semialgebraic map whose image is the semial- 
gebraic arc C. Set A := Lp{a) and B := ip{h) so that C connects A to 
An ordered sampling of C is an ordered collection of points 

y := {Pi,...,P„} C c, 

such that the collection 

:= {ti=^~\Pi),...,tn ■.= ^~\Pn)] C [a, 6] 

satisfies 

ti <t2 < ■■■ <tn. 

The mei/j of the ordered sampling T is the positive number 

||y|| := max{dist(A, Pi), dist(Pi, P2), . . . , dist(P„_i, P„), dist(P„, B) }. 

We denote by C{'J') the PL-curve (Pi, . . . , P„). We will need the following result whose proof is 
delegated to Appendix A. 
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Proposition 2.8. Suppose that C is a semialgebraic arc and for very e > we are given and ordered 
sampling 7^ ofC. Denote by L (resp. K) the length (resp. total curvature) ofC and by (resp. Kf,} 
the length (resp. total curvature) of Oil's)- If 

lim liyJI = 0, 



then 

£\0 " e\0 



lim Ls = L and lim = K. 



□ 



Theorem 2.9. Suppose that h : [a,b] — )■ K is a continuous semialgebraic function and e — )■ cr(e) is a 
spread function satisfying the additional condition 

lim ea{ef = oo. (2.4) 

For every e > we choose an e-sample with spread a"(e) of the graph T of h. Denote by the 
graph of the PL-function L^^ described in Definition 2.3. Then, as e \jOwe have 

length (Ce) ^ length (r) and K{Ce) ^ K{T). 



Proof. We use a simple strategy. More precisely for every e > we construct an ordered sampling 
of r such that 

lim||yj|=0, (2.5) 
and if denotes the PL-curve T{'S's) determined by the ordered sampling Ve, then as e \ we have 

length (C) = length (r,) + O , (2.6a) 

K{Ce) = K{Te) + O (^^) ■ (2.6b) 

The desired conclusions will then follow from Proposition 2.8. 

Suppose that consists of the points Qq, Qf , . . . , Q^j-^) arranged in the increasing order defined by 
their x-coordinates. Observe that since has spread a{s) then 

nie) < ^ (2.7) 

Each of the points of is the center of a pixel that touches P. Thus, for any A; = 0, 1, ... , n(e) there 
exists a point G P that lies in the same pixel as Q^. We obtain in this fashion an ordered sampling 

= {^0 1 A^' • • • ' Pn(e) } 

of P. The function h is continuous and semialgebraic and thus it is Holder continuous with some Holder 
exponent a £ (0, 1]. This proves that 

m\=0{{ea{e)r). 

The condition (2.5) now follows from the property (2.1b) of a spread function. 
From the choice of the points we deduce that for any A; = 1, . . . , n(e) we have 

-eV2 < dist(Pfc^_i, P|) - dist(g|„i, Ql) < eV2 
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SO that by summing over k we deduce 

^^^-^ — ^ < -n(e)eV2 < length (rj - length (C^) < n(e)eV2 < \ , — ^. 

cj(e) cr(e) 

The equahty (2.6a) now follows from the property (2.1a) of a spread function. 

Now we turn to total curvature. For k = 1, . . . ,n{e) we denote by ?7i|. the slope of the segment 
[P|„j^, P|] and by m| the slope of the segment Note that for any e > and any k = 

1, . . . , n(e) we have from the definition of a spread that 

dist(Q|_i,Q|) <e<T(e). 



Furthermore we have shown that 
These two inequalities imply that 



distiQ%,Pl)<eV2. 

\m% — m|| = O 
There exist 6*1, ^| G (-f , f ) such that 



m% = tan 0| , m| 
The difference formula of tangent impUes that: 



9| — 0| = arctan 



1 + m|m|. 



Using the above equation and the fact that |m| — m|| = O y^^j > we see that there is a constant C 
independent of e such that 

C 

' " a{e) 

and therefore 

' ' cr(e) £a[eY 

The equality (2.6b) now follows from (2.4). 

□ 

Corollary 2.10. Let S{l3, r) be an elementary set. Fix a spread a{e) function satisfying the condition 
(2.4). For each e we choose compatible e-upper/lower profiles r.f- of S with spread cr(e). We denote 
by Se the PL approximation of S defined by these samples. Then 

lim length {dSe) = length {dSe) (2.8a) 

e\0 

\\mK{dSe) = K{dS). (2.8b) 

Proof. The semialgebraic functions /3 and r are differentiable outside a finite subset of (a, b) and the 
limits 

/3'(a) := lim /3'(x), r'(a) := lim r'(x), 
/3'(6) := lim/?'(x), t'(6) := lim r'(x) 



exist in [— oo, cx)]. Let 
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The compatibility condition implies that 

^k=^k--^k, < ' V/c = 0, 1,... ,n. 

Let /3j be the PL function whose graph is L^- and be the PL function whose graph is L.^+. Let 

mf (e) indicate the slope of the z-th line segment of the graph of (3^ and similarly let mj{e) indicate the 
slope of the i-th line segment of the graph of Tg. We have 

length(55,) = length (F^J + length (F^J + dist(eo", C^) + dist(q,), )• 

Theorem 2.9 implies that as e \ we have 

length (r^jj ^ length (F^), length (T^J ^ length (T^). 

Moreover, as e \ we have 

dist(^o-, ej) ^ dist(/3(o), r(a) ), dist(q^), ) ^ dist( /3(6), t(6) ). 

This proves (2.8a). 
Similarly 

K{dSs) = Itt — arctan(m^(e) )| 

n 

+ I arctan(m^(e) ) — arctan(m^_j^(e) )| 

j=2 

+ I arctan(m^(e)) — vrj + jvr — arctan(my)| 

n 

+ I arctan(77T,J(e) ) — arctan(mj„]^(e) )| 

i=2 

+ I arctan(m^(e) ) — 7r| 

which can be rewritten as 

K{dSs) = Ivr — arctan(m^(e) )| + | arctan(m^(e) ) — vrj 

+ Itt — arctan(my(e) )| + | arctan(m^(e) ) — tt\ (2.9) 

Theorem 2.9 implies 

lim K{(3,) = K{(3) and lim K{Te) = K{t). (2.10) 

Now note that each line segment is defined by connecting two points in the pixelation of /3 or r over an 
interval of at most £a{£). As e — )• we have 



lim (e) = p' ia) , lim ie) = /?'(&), 
lim?Ti[(e) = T'(a), limm^(e) = /3'(6). 

e\0 e\0 



(2.11) 



Combining (2.9), (2.10), (2.11) we find that 

lim K{dSe) = Ivr — arctan(/3'(a))| + | arctan(/3'(6)) — 7r| 

+ Ivr — arctan(T'(o))| + | arctan(T'(6)) — 7r| (2-12) 
+ K(/3) + K(r). 

Note that [vr — arctan(/3'(a))| is the value of the angle between the vertical line x = a and the tangent 
line to the graph of /3 at (a, /3(a)). Similarly each other difference on the right hand side of (2.12) 
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corresponds to an angle at one of the corners of dS. Therefore the right hand side of the (2.12) is equal 
to the K{dS), so the corollary holds. □ 

3. Separation results 

In the previous section we have dealt only with the elementary regions and we investigated mainly 
geometric properties of these regions and their pixelation. In this section we turn our attention to the 
relationship between the topologies of a semialgebraic set and those of its pixelations. 

Surprisingly, this is a nontrivial matter In general, the homotopy type of a planar set may be quite 
different from those of its pixelations. This can happen even for a simple PL case. Consider the set 
S composed of the rays starting from the origin and proceeding in the positive direction with slopes 
^ and ^ (see Figure 6). A careful examination of P£{S) reveals the existence of cycle. Worse yet, 
the pixelations for smaller values of e are simply rescalings of the larger pixelations. This means that 
PeiS) contains a cycle for all small e. For a taste of how much worse this can get we refer to [12, 
Appendix A]. 



Figure 6. The pixelation of the lines with slope ^ and After this point the pixela- 
tions of the upper and lower lines diverge permanently and no more holes are formed. 

The next result provides a first ray of hope. For any compact set X C we denote by S(X) the 
set of connected components of X. 

Proposition 3.1. Let S cM? be a compact semialgebraic set. Then for sufficiently small e, the number 
of connected components of Pe{S) agrees with the number of connected components of S. 

Proof. We have a natural map 6(5) — )• C( Pe{S) ) that associates to each connected component C of 5 
the unique connected component of PeiS) containing C. For e sufficiently small this map is injective. 
Since Pe{S) contains only pixels that intersect S, we deduce that each connected component of Pe{S) 
contains at least one connected component of S. Thus, PeiS) has at most as many components as S. □ 

The above result guarantees that the zeroth Betti number of a compact semialgebraic set coincides 
with those of its sufficiently fine pixelations. Proposition 3.1 also suggests that, for small s, the only 
way that the homotopy type of -Pe(S') can disagree with that of S is if Pe{S) has holes, i.e., cycles in 
PeiS) that are not contained in the image of the inclusion induced morphism 

Hi{S,Z)^Hi{P,{S),Z). 

Thus recovery of S from Pe{S) will depend on distinguishing the cycles of Pe{S) that correspond to 
real cycles from S from those that are merely artifacts of the pixelation. To discard these holes, we 
adopt a strategy inspired from Morse theory. 
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Definition 3.2. Let 5 C be a compact set and e > 

(1) For every e-generic x we set 

n£{x) = ns,e '■= # of connected components of CelS, x). 

(When the set S is understood from context we use the simpler notation instead of ns^s-) 
We will refer to ns{x) as the stack counter function of S. 

(2) For any xq G M we define 

n{xQ) = ns{x) := # of components of the intersection of S with the vertical line x = xq- 

We will refer to ng the component counter of S. 

(3) K jumping point of ns is a real number xq such that 

ns{xQ) / nsixQ) := lim ns{x) or ns{x) / ns{x^) := lim ns{x). 

X/^Xq x\X0 

We denote by 3 s the set of jumping points of ns- We will refer to 3 s as the jumping set of S. 

(4) A jumping point of is ai^eal number xq G + f such that 

ne{xo - e) / n(xo). 

We denote by 3s,e the set of jumping points of £. We will refer to it as the e-jumping set of 
S. 

□ 



Let us point out that if S is semialgebraic, then its jumping set 3s is finite and it is contained in the 
set of critical values of the restriction to S of the function h{x, y) = x. The function tells us how 
many stacks are in a column. The jumps of are a first indicator of the presence of cycles in Pe(5'). 
To decide whether they ai^e holes, as opposed to cycles coming from S we will rely on the next key 
technical result. 

Theorem 3.3 (Separation Theorem). Let f,g : [a, b] M be two semialgebraic continuous functions 
such that f{x) < g(x), Vx E [a,b]. Denote by G the union of the graphs of f and g. Fix L > 0, 
a S (0, 1] and xq S [a, b] such that either 

\g{x)-g{y)\<L\x-y\'', or \f{x) - f{y) < L\x - y\''yx,y € [a,b] (3.1) 

and 

9{xo) - f{xo) < g{x) - f{x), VxG[a,6], (3.2) 

Then for any e > such that 

3e + Le" < ^(xo) - /(xo) (3.3) 
and any e-generic x G [a, 6] \ eZ the column Cs{G, xq) has two components. In other words, if 

min ( g{x) — /(x) ) > 3e + Le", 

x(^[a,b] 

then for any e-generic x € [o, b] we have 

nG,e{x) = ncix). 



Proof. We deal with the case that \g{x) — g{y)\ < L\x — Note that 

Cs{G,x)=Cs{f,x)uCe{g,x), 
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and furthermore, Theorem 1.3 impHes that each of these columns is connected. Therefore Ce(G, x) 
will have two components exactly when Ce{f, x) and Ce((7, x) do not intersect. Since f < g and x is 
e-generic, this will occur when 

T,{f,x)<B,{g,x) 

or equivalently, 

B,{g,x)-T,{f,x)> e. 

Now fix j; e [a, b] and let i G Z such that ie < x < {i + l)e. Choose Xf,Xg e [ie, {i + l)e)] such that 
f{xf) = max /(x), g{xg) = min g{x). 

Therefore we have 

Be{g,x)>g{xg)-e, Te{f,x) <f{xf) + s, 



so that 
Thus 



Be{g,x)-Te{f,x)>g{xg)-f{xf)-2e. 

Be{g, x) - Ts{f, x) = g{xg) - g{xf) + g{xf) - f{xf) - 2e 
(3.2) 

> 9{xg) - g{xf) + g{xo) - f{xo) -2e 

(3.1) (3.3) 

> g{xo)-f{xo)-Le''-2e > e. 

which completes the proof for the case that \g{x) — g{y)\ < L\x — ?/|". If instead we are given that 
\f{x) — f{y)\ < L|x — the proof is similai^, except that the roles of / and g exchanged. □ 

This important Separation Theorem can be used to prove the following two results, which will tell us 
exactly when n and correspond, using only information from the pixelation. With these theorems 
we will be able to distinguish real cycles from the original set from fake cycles created by the pixelation. 

Theorem 3.4. Let S C M? be a compact semialgebraic set with jumping set 3^. Then there exist 
1^0 = ^o(5') G (0,1], 1^0 = J^oiS) > and Eq = eo{S) > 0, depending only on S, such that, if 
< e < eo X is an e-generic value such that 

dist{x,3s) > i^oe''", 

then ns^e{x) = ns{x). 

Proof. Let xq < xi < • • • < be the jumping points of n = ng. We set 

Axj := Xj — Xj-i, Vi = 1, . . . A := min Axj. 

l<i<i 

Note that n(x) is constant on each of the intervals (xi_i, Xj). For i = 1, . . . , ^ we set 

Si := {{x,y) £ S; x G [xj_i, Xj] }. 
The set Si is a disjoint union of elementary sets (see Definition 2.5 for notations) 

S{l3i,j,Tij), j = 0, . . . ,pi, 

"stacked one above the other", i.e., 

/3i,o{x) < Tifi{x) < Pi,lix) < Ti,l(x) < • • • < /3i,p^(x) < Tj,p,(x), Vx e (Xi_i,Xi). (3.4) 

From Proposition 2.6 we deduce that for any e-generic x € (xi_i, Xj) we have n(x) = pi. 

Both of the functions /3j j and Tij are continuous and semialgebraic. For any i = !,...,£, any 
j = 1, . . .pi and any h G (0, ^) we denote by jij{h) the minimum of f3ij — nj-i on the interval 
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[xi-i + h,Xi — h]. From Lojasewicz's inequality we deduce that there exists C = C{S) > and 
r = r{S) G Z>o such that for any i = 1, . . . ,£, any j = 1, . . . ,pi and any h £ (0, ^) we have 

7,,,(/i) > Cff. (3.5) 

Fix L > and a > such that for any i = I, . . . ,£, j = 1, . . . ,pi and any x, y G [xi^i, Xi] we have 

- + - rij{y)\ < L\x - y\°'. 

Now choose Eq > 0, uq > and kq G (0, 1] such that 

C{uoe^''Y >3e + Le'', VeG(0,eo)- (3-6) 

The conclusion now follows from (3.5), (3.6) and Theorem 3.3. □ 

Definition 3.5. The constant kq{S) guaranteed by Theorem 3.4 is called the separation constant of the 
set 5. 

This theorem tells us that the jumps of occur within uqe'^^ pixels from the jumps in n. A priori, 
it could be possible that, given a jumping point xq of n, there is no jump in within uqe'^" pixels of 
xq. This next theorem shows that in fact this cannot happen. 

Theorem 3.6. Let S be a generic compact semialgebraic set, and Eq = £oiS), = i^q{S) as in 
Theorem 3.4. Let kq = i^oiS) be the separation constant of S. Then, there exist £i = £i{S) > such 
that if £ < min(eo, £i) and xq is a jumping point ofn = ns, then = ns^e h'^^ <^t least one jumping 
point in the interval [xq — i'o£'^° ■, xq + foe'^"]. 

Proof. Fix a good stratification 3" of 5 such that the function h{x, y) = x is a stratified Morse function 
with respect to {S, 3"). Then there exists exactly one critical point po G 5 of /i such that h{pQ) = xq. 

Letpo = (a;o,yo)- 

Since xq is a jumping point of n we have 

n(X(J") / n(xo) or n(xo) / n(xg ). 

We discuss only the case n(xQ ) / n(xo) because the other case reduces to this case applied to the 
region obtained from S via a reflection in the y-axis. For every e > fix an e-generic point x'q{£) such 
that 

xo, if xo G M \ 
Xq — I if Xo G eZ. 
We distinguish several cases cases. 

Case 1. n{x^) > n(xo). We can find 5 > sufficiently small such that the interval [xq — 5, xq] will 
contain no new jumping points of S. The set 

S[xo~SM ■■= {{x,y) £ S; X G [xo - 5,xq\ } 

is disjoint union of elementary regions 

S{/3j,Tj), j = 0, . . . ,m = n(xo ) - 1, 

"stacked one above the other", i.e., 

/3o{x) < To(x) < /3i(x) < Ti(x) < • • • < /3„,(x) < Tm{x), Vx G (xq - 6,Xo), 

where Pj,Tj are are continuous semialgebraic functions. Since n(xo) < ■".(xq ) we deduce that there 
exist jo, ji = 1, . . . , m such that jo < ji and 

/3ji(xo) = Tjo-i(xo) and := /3j(xo) - rj_i(xo) > 0, Vj [jo, ji]. 



Xo(e) 
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In particulai^, for any e > 0, the e-stacks over XQ{e) of 

^Wji^'^ji)^ • • • ^ S{j3jQ,TjQ)^ S{l3jQ-i,Tj^-i) 

have a point in common. 

Now choose ei sufficiently small so that for j ^ [jOiJi] and e < ei, the e-stack of S{/3j,Tj) over 
Xq^e) is disjoint form the e-stack of S'(/3j_i, rj_i) over XQ(e). Fix e < min(eo,ei). The above 
discussion shows that 

n(xo) = neix'oie)). 

Theorem 3.4 now implies that 

ns{xo - foe''°) = n{xo - i^oe"") > n{xo). 
This proves that the interval [xq — uqe'^^^xq] contains a jumping point of n^. 

Case 2. n(xQ ) < n(xo). The critical point po has the property that there exists a tiny disk D centered 
at pq such that the intersection of D with the open half -plane < xq} C is empty. In particular, 
this shows that po is an isolated point of the set 

Sx<xo = {ix,y) ^ S; x<xo}- 

If n{xQ) = 0, the conclusion is obvious. We assume that n(a;Q ) > 0. Choose 6 > such that the 
interval [xq — S,xo) contains no jumping point of S. Set 

R := cZ(5[^^_5_^.q] \ {po}). 

Then i? is a union of simple regions 

S{l3j,Tj), j = 0,1..., m = n{x^) - I, 
where Pj and tj are continuous semialgebraic functions functions such that 

Po{x) < To(x) < (3i{x) < ri(x) < • • • < /3m{x) < Tmix), Vx G [xq - S,Xo]. 
We can find ei = ei {S) such that for any e < ei and any e-generic x G [xq — S, xq] we have: 

• UR^eix) = ns{x) = m + l = ns{x^), and 

• the e-column of S^xq-s^xq] ov^i" consists of n(xo) = m + 2 stacks. 
Theorem 3.4 implies that 

ns^eix) = ns{x) = ns(xQ ) = m + 1 Vx G [xq - (5, xq — z^oe""] \ ^e- 

On the other hand, ns^ei^Q) = m + 2. Thus the interval [xq — uqe'^^ , xq] must contain a jumping point 
of ns,e- 

□ 

Remark 3.7. Theorem 3.4 states that the two functions n and coincide at points situated at a 
distance at least i'Q{S)e'^°~^ pixels away from the jumping points of n. On the other hand, Theorem 
3.6 shows that, for a generic semialgebraic set, then within ijQ{S)e'^"~^ pixels from a jumping point of 
n there must be jumping points of n^. □ 

Definition 3.8. Let S be a generic semialgebraic set in and the constants eq (S) and ei (5) as defined 
in Theorems 3.4 and 3.6. We set 

h{S) ■.= mm{eo{S), ei{S)), 
and we will refer to it as the critical resolution of S. 

□ 
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4. Approximation of generic semi-algebraic sets 

This section is the heart of the paper. Here we will create an algorithm which will approximate a 
generic semi-algebraic set using only its pixelations, and then prove a very strong convergence result 
for this approximation. This algorithm is based on the central algorithm of [12], updated to handle the 
additional complexities of semi-algebraic sets. 

We first observe that when narrow vertical strips around the jumping set Js are removed from a 
semi-algebraic set S, the remainder is a disjoint union of elementary sets. Corollary 2.10 indicates a 
good way to approximate continuous semi-algebraic functions, and Theorems 3.4 and 3.6 indicate that 
for small e, the jumping points of S become close to the jumping points of Pe(5). Therefore a viable 
approximation technique is to treat parts of Pe{S) which occur near jumping points as noise (to be 
approximated crudely) and to approximate outside of this noise by means of Corollary 2.10. 

There are two quantities which must be used in this approximation. The first is the previously 
mentioned spread function a which determines the width of line segments to be used in approximating 
outside of noise. From Corollary 2. 10 we know that this spread function should satisfy the following 
limits: 

lim ea{e) = 0, lim e(cr(e))^ = oo. 

e\0 e\0 

The second quantity determines the width of the noise, measured in pixels, about jumping points. We 
will call this quantity u and refer to it as the noise width. It is defined as follows: 

Definition 4.1. Let 5 be a semi-algebraic set and kq be its separation constant. Then a noise width v 
is a function v : — )• Z"*" which satisfies the following equations: 

limei/(e)=0, (4.1a) 

lim = oo, (4.1b) 

The first property in this definition ensures that noise is a small phenomenon. The second property 
implies that eu{£) increases faster than e'^° so that the noise will eventually contain all fake cycles 
(consult Remai^k 3.7). 

For a reasonable approximation we must have an estimate of kq. For example if kq = ^ we could 
set v{e) = [e~3"]. 

Using these constants together with the techniques developed in the first three sections we create an 
algorithm to approximate a semi-algebraic set from its e-pixelations. For ease of discussion we will 
use the following language in the algorithm: 

Definition 4.2. Let 5 be a semi-algebraic set and Pe{S) its e-pixelation. If A C M then the part of 
over A is the set 

Pe{S) n (A X M) 



Algoritlim 4.3. (1) Choose a spread a such that ea{eY — ^ oo and e(T{e) — )■ as e — )• 0. 

(2) Choose a noise width u such that — > oo and ev{e) — )■ as e — )• 0. 

(3) (a) For each point p ^ let C~ (p) be the largest odd multiple of | such that p — C7 (p) > 

ev{e). Similarly let Ctip) smallest odd multiple of | such that p + Ctip) < ^^^C^)- 

(b) For each p & Je let Ae(p) be the interval [(~{p),C^{p)]. 

(c) Let 

A, = U A,(p). 
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The set Ag is called the noise set of Pe{S) and its connected components are called the 
noise intervals of Pe{S). 

(4) Define IR^ to be the closure of M \ A^. We call 'Ji^ the regular set of Pe(5), and its connected 
components are called regular intervals of Pe{S). 

(5) For each bounded regular interval / C and each connected component C of Pe {S) H (/ x M), 
the part of Pe{S) over the regular regular intervals / do the following: 

(a) Choose compatible upper and lower profiles 11+ and IT" on the component C. 

(b) Choose compatible upper and lower samples and H~ with spread a. 

(c) Generate the PL-approximation determined by the above upper and lower samples. 

(6) The union of all PL-approximations found in the above step is called the regular approximation 
which we denote by 5/"^^. 

(7) For each noise interval / C A^ denote by Ce(/) the set of connected components of Pe{S) 
over / 

(a) Let / = [a, h] C A^ be a noise interval of Pe{S). 

(b) For every C G Ce{I) we denote by Uc (resp. Lc) the highest (resp. lowest) y-coordinate 
of the center of a pixel in C. 

(c) Denote by IPe(C) the rectangle determined by the inequalities 

X > a,x <b,y > Lq, y <Uc 

(see Figure 7). 
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Figure 7 . Covering up the noise. 



(d) The noise approximation over I, which we denote by Pe{I), is the union 

PAI):= U ^eiC) 

C&Ce{I) 



(8) The union of all Ni;{I) where / are the noise intervals in Ag is called the noise approximation 
which we denote by S^°^^^. 

(9) The final approximation is simply the union of the noise and regular approximations, i.e. 



g g noise 1^ g reg 

This final set Sg will be piecewise linear by construction, and will be a good approximation of the 
set. □ 



Here "good approximation" means that it captures both topological and geometric information such 
as area, perimeter, and curvature measures. The precise notion of "good approximation" relies on the 
concept of normal cycle. 
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The normal cycle is a coiTcspondence that associates to each compact planai" semialgebraic set X a 
1-dimensional current N-^ on the unit sphere tangent bundle of 'M? 

For a precise definition of this object we refer to [1, 6, 9, 10]. Here we will content ourselves with a 
brief description of its construction. 

For a semialgebraic domain D C with -boundary the normal cycle has a simple descrip- 
tion. It is the current of integration given by the closed curve C S{TB?) 

gz) = {(n(p),p)G5(rM2); p^dD] 

where n(p) denotes the unit outer normal to dD at p G dD. Equivalently, Sd is the graph the the 
Gauss map 

dD Bp^n{p) G S^. 

Clearly, in this case, the normal cycle contains all the curvature information concerning the boundary 
of D. 

More generally, if 5 is a compact semialgebraic set, the we can find a semialgebraic function 
/ : ^ [0, oo) such that S = f~^{0). For all e > sufficiently small the region Se := {/ < e} 
is a semialgebraic domain with C^-boundary so we can define the normal cycle N^^ as above. One 
can show that as e — )• the currents N^^ converge weakly to a cuiTcnt which by definition is the 
normal cycle of S. The hard part is to prove that this current is independent of the choice of defining 
function /. This cuiTcnt is a cun^ent of integration along a finite number of oriented semialgebraic arcs 
in S'(rM^). We refer to ["-J] to a more in depth description of the normal cycle of planar semialgebraic 
sets. In particular, in [9] one can see how this current captures the various curvature properties of S. 

Theorem 4.4. Let S be a generic compact semi-algebraic subset of the plane. For each e, let be the 
PL set constructed through the above approximation algorithm. Then converges to N^^ weakly, 
where and N^^ indicate the normal cycles of S and respectively. 

Proof. First we note that the approximation converges in the Hausdorff metric to the original set. This 
is because each vertex of a line segment is taken from a pixel which contains a piece of the boundary 
of the original set. Since every pixel of the pixelation can be at most e\/2 far from the original set, this 
forces the approximation into a tube around the original set which becomes arbitrarily small as e goes 
to 0. 

From here the strategy of the proof will make heavy use of the inclusion-exclusion principle satisfied 
by the normal cycle correspondence X i— )• . More explicitly, this means that for any compact semi- 
algebraic sets X and Y, we have 

We will use this principle to reduce the calculation of the normal cycle of into calculations on 
simpler subsets of S. 

First, we need to introduce some more notation. We set 5*0 = S. For each e > and each c G J s,e 
we indicate by Je(c) the e-noise interval containing c. For each c G J s we set Jo(c) := {c}. We then 
define, for each e > 0, the noise strip J^eic) as 

N,(c) :={(x,y);x G Je(c)} 

finally we set 

and -Re as \ ^ST^. 
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For each e > we construct a graph as follows. The vertex set is the set of connected components 
of T^^nSe. The edge set is the set of connected components of OlsCiSe, so that two vertices are connected 
if and only if there is a component of n Se which connects the two appropriate vertices of 3sfe n S'e. 
This graph is the Reeb graph of the projection of Se onto the .x-axis. Furthermore there is a 5 > 
independent of S such that for all e G (0, 6] the graph is isomorphic to the graph Fq. Let 

£2 := Tam{h{S), 5} 

where h{S) is the critical resolution defined in Definition 3.8. For the remainder of the proof we will 
deal only with e G (0, £2]- 

Let be the set of vertices of F^ and E.^ be the set of edges of Fg. Note that, by the above 
equivalence of Reeb graphs, for e G [0, 82], there is a natural bijection between the vertices of F^ with 
those of Fq and similarly for the edges. For any vertex v of Vq and e G [0, £2] we indicate by ^ 
the connected component of n Ss corresponding to the vertex. Similarly, for any edge e of Fq we 
indicate by Ce,e the closure of the connected component of Ol^ n corresponding to e. We have the 
following result, (compare [12, Lemma 5.3]). 

Lemma 4.5. For any e G [0, £2] H'e have 

jY^e _ jyC^.e _|_ ]\fCe,e _ ]\^C^^enCe,e _ (4.2) 

vGVq eG£o veVo ee£() 

Proof. Note that we have a decomposition 

5, = ( I J , 1 U I \\Ce,s\. (4.3) 




We need to discuss separately the cases e > and e = 0. 

1. Assume that e G (0, £2]- In this case we have 

C.,.nC„.,, = = Ce,enCe',e, yv^v', 6 / e'. (4.4) 

The equality (4.2) now follows from inclusion-exclusion principle applied to the decomposition (4.3) 
satisfying the overlap conditions (4.4). 

2. Assume that e = 0. In this case the overlap conditions are more complicated. We have 

Cv,o n C„/^o = 0, V?; / v', (4.5a) 

Ce,o n Ce',0 = ^ e n e' = 0, (4.5b) 
where the condition e n e' = signifies that the edges e and e' have no vertex in common. Moreover, 

fl Ce,o = C,,o, G Vo, Ac E,. (4.6) 

eGA 

Using (4.3), (4.5a), (4.5b), (4.6) and the inclusion-exclusion principle we deduce 

jV^ _ ^ jsfCvfi _|_ ^ jyCe.o _ ^ ^ jyC^.onCe.o 
i>GVo eG£o veVo eeE^ 

+ (_l)l^l + lAr^-n{ag^)Ce,0 + ^ (_l)l^lArn.eAC^e,0 
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veVo ee£o uGVo ee-E„ 

+ ^ ( y: ( (-1)1^1^^+ (-1)1^1 )) N^-' 



□ 



The above lemma means that the proof will follow if the following three equations are satisfied for 
every edge e and vertex e in : 

lim N^^-' = Af^"'" (4.7a) 
lim JV^''.£^'^e,e _ ]yC^,o<~^Cefi (4.7b) 

lim N'^''^'' = N^'^'O (4.7c) 

where convergence of each limit is meant in the weak sense of currents. 

Each of these equations will rely on an approximation result of normal cycles proved by Joseph Fu 
in [^]. A restricted version of this theorem, which shall suffice for the purposes of this paper, is stated 
below: 

Theorem 4.6 (Approximation Theorem). Suppose S is a compact semialgebraic subset of the plane 
and for each e > we are given a compact semialgebraic subset of the plane with the following 
properties. 

(1) There is a compact set K C which contains each S^- 

(2) There is a M G M such that 

mass( A/"^^ ) < M, Ve. 

(3) For almost every ^ G Hom(]R^, M) and almost every c G R we have 

\\m x{Se n > c}) = x{S n > c}) 

Then converges to N as e ^ weakly and in the flat metric. □ 

Equations (4.7a) and (4.7b) will both follow from applying this approximation theorem to the case 
of rectangles. Therefore, the following lemma will be useful: 

Lemma 4.7. Suppose {Ss)e>o a family of convex polygons in the plane that converge in the Hausdorff 
metric to a convex polygon S. Then N^^ converges weakly to as e ^ 0. 

Proof. We argue by proving the conditions of Fu's Theorem. Observe first that there exists R > such 
that 

dist{Se,S) < R, Ve 

and thus the condition (1) of the Approximation Theorem. The computations of P, Chap. 23] show 
that mass of the normal cycle of a convex polygon P is equal to 27r + L{P). From Hadwiger's charac- 
terization theorem [ , , Thm. 9.1.1] we deduce that 



limL(5e) = L{S) 



and thus condition (2) is also satisfied. 
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Therefore we must show that for almost every ^ G Hom(]R^, M) and almost every c G M we have 

lim y(S, n < c}) = x{S D < c}) 

Note that 5 and each Se ai^e all convex subsets of the plane. Therefore any intersection with a half -plane 
will either be empty or be a contractible set. Therefore to prove the convergence of Euler characteristic 
on half-planes we need only prove that a half plane H will only intersect Se for small e if and only if it 
intersects S. This is true since H Ci Se converges in the Hausdorff metric to H n S. □ 

Proof of (4.7a) Fix a vertex v G Vq. The set C^fi is a subset of a vertical line over a jumping point, 
and is so either a point or a line segment. For every e G [0, £2], the set C^^e is a rectangle which spans a 
noise interval and contains C^^o- The width of (and so the width of C„,e) is proportional to ez^(e), 
and so vanishes as e — (by choice of v). 

The rectangle Cv,e is constructed by choosing the highest and lowest pixels from the component of 
le n Pe{S) containing Cv,o- For sufficiently small e, the noise interval le will be thin enough so that 
Cv,e n S can be described as a number of regions lying between the graphs of functions which are 
everywhere except possibly at the jumping point. This implies that for small e, the height of ^ 
differs from the height of Cefi be an arbitrarily small amount. 

Since the height and width of Cv,e converge to the height and width of C^fi and since each ^ 
contains C^,o> the rectangles C^t, £ converge in the Hausdorff metric to C^,o- Therefore by Lemma 4.7 

Proof of (4.7b) Fix a vertex e G £0 and a vertex v G Vq. If Ce,o H C^fl = 0> then for sufficiently small 
e the component Ce,e will also not intersect Cv,e and so the convergence in normal cycle follows. 

If the Ce,o and C^fi do in fact intersect, then note that Ce,o H C^fl = C^fl (since the vertex is a 
connected component over a point, and the edges is a connected component over an interval which 
overlaps that point). 

The intersection Ce,e H C^.e is a vertical line segment. In fact it is either the right or left edge of 
C-i,,e- However, since C^^e converges to a vertical line segment Ce,o> it follows that its left right and 
right edges converge to the same line segment. Therefore (4.7b) follows from (4.7a). 

Proof of (4.7c) We again plan to use the Approximation Theorem. Condition (1) of the theorem is 
plainly satisfied while condition (2) follows from Corollary 2.10 and the explicit description of the 
mass of the normal cycle of a planar set given in [ ■, Chap. 23]. All that is left to do is to verify 
condition (3) of the Approximation Theorem. 

The component Ce,o is an elementary region defined by continuous semialgebraic functions 

More precisely, this means that 

/3e{x) < Te{x), Vx G [o,6], 

and 

Ce,o = G x£[a,b], f3e{x) < y < Te{x) ] . 

There exists an integer ?i > and points 

a = Co < ci < ■ ■ ■ < Cn = b 

such that for any i = 1, . . . , n the restrictions of /3e and Te to (cj_i, Cj) are real analytic. Moreover, 
since the set S is generic, the derivatives /3g and are bounded near ci, . . . , c„_i. In particular, the 
functions (3^ and are locally Lipschitz on the open interval (a, b). We will refer to the points 

(cj,/3e(Cj)), {cj,Te{cj)), j = 0, l,...,n, 
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as the vertices of Ce,o- The other points on these graphs ai^e called regular. Now fix a constant c G R 
and a linear map ^ : — )• R, y) = ux + vy, with the following generic properties: 

Gi. The restriction of to Ce,o is a stratified Morse function, and f / 0, i.e., the level sets of ^ are 

not vertical lines. 
G2. The constant c is not a critical value of ^Ic^ q. 

G3. The line c := {£, = c} does not contain any of the vertices of Ce,o- 
We will show that 

limx(Ce,en{e> c}) =x(t^e,on{e> c}). (4.8) 

We denote by C| q the intersection of Ce,o with the vertical strip 

Oe := a + iy{e)e < x < b — u{e)e =: bg. 

for all e > small enough so that the above inequalities make sense. Let us observe that the conditions 
Gi,G2 and G3 imply that for e sufficiently small we have 

X ( a,o n {e > c } ) = X ( C|,o n > c } ) . 

Thus, to prove (4.8) it suffices to show that 

x(Ce,.n{e>c}) =x(c|,onU>c}), Ve«l. (4.9) 

The region Ce,e is an elementary region defined by the PL functions 

Pe,e,Te,e ■ [ae,be] M, f3e,e{x) < Te,e{x), Vx G [ae,be]. 

Here we develop some terminology to handle the intersections of these PL boundary functions. If 
/ : [s, t] — )• M is a piecewise function, then we say that the line L^ c intersects the graph of / 
transversally at a point Pq = (xo, f{xo)) if there exists a.6 > such that the function 

is differentiable on the set < |x — xo| < 6 and it has constant sign on this set. We will denote by 
sign(pQ, /) G {±1} this sign. Thus, if sign(po, /) = 1, then the curve 

X I—)- (x, /(x)), \x — xol < 6, 

intersects the line L^ c at Pq coming from the half-plane < c} and entering the half -plane > c}. 
For any point Pq = (xq, yo) and any r > we denote by T,r{p) the closed square 

Sr(p) := {{x,y) G M^; |x - Xol, |y - yol <r}. 
We need to discuss separately three cases. 

Case 1. The elementary set Ce,o nondegenerate, i.e., fie{x) < Te(x), Vx G {a,b). The conditions 
Gi,G2,G3 imply that there exists a compact subinterval / = [a*, 6*] C (a, 6) such that the the 
line c intersects the graphs of /3e and Te transversally in regular points on these graphs whose x- 
coordinates are contained in the interval /. Denote by /J3 (resp. J^) the intersection of L^ c with the 
graph of Pe (resp. Tg.) The superscript of in /° comes our convention that S = Sq. Finally we set 

Fix a small positive real number r with the following properties (see Figure 8). 

• The closed squares T,r{p), p G -T*^ are pairwise disjoint. 

• For each point p G there exists i = 0,1, . . . ,n such that the projection of Sr(p) onto the 
X-axis is contained in a compact sub-interval Jp C (cj_i, Ci). 
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Figure 8. Isolating the intersection points of L^^c with the graptis of j3 and t so that 
the squares Sr(p) do not touch any of the vertical lines containing the singular points 
of these graphs. 

^Vc set 

Lemma 4.8. (a) Denote by I^^ the intersection of L^ ^. with the graph of /3e,e- There exists > with 
the following properties. 

(al) For any e < epwe have 

(a2) For any p ^ and any e < £p the line L^^c intersects the portion of the graph of j3e^£ inside 
Dr{p) in a unique point p(e). This intersection is transversal and 

sign(p,/3e,o) = sign(p(e),/3e,e). (4.10) 
(b) Similar statements are true with the bottom functions (3^^^ replaced with the top functions Te^e- 

□ 

We defer the proof of this result to the end of this section. 

Set = U We will refer to the intersection of L^ ^. with 9C| q as the 0-crossing set and, 
for e > 0, we we will refer to the intersection of L^ ^. with dCe,e as the e-crossing set. For e > we 
denote by the e-crossing set. 

Observe that the sets is contained in but the e-crossing set may contain additional points, 
namely, the intersection of L^ ^. with the vertical lines x = a^,be. The Hausdorff distance between the 

and \q goes to zero as e \ 0. Moreover, Lemma 4.8 implies that for any e sufficiently small there 
exists a bijection 

: to ^ t. 

defined by 

= t,ns,(p). 

For e > we denote by C^^ the intersection of Ce,e with the half -plane > c}. Similarly, we define 
Cq to be the intersection of C| q with the same half-plane. We have to prove that 

X(C+) = x(Co+), Ve«l. (4.11) 

For e > the connected components of C+g are all homeomorphic to closed 2-dimensional disks so 
that the Euler char acteristic of is equal to the number of boundary components of dC^. 
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Observe that if \q = 0, then fe = for all e sufficiently small. In this case C^f^ is homeomorphic 
to a closed disk for all e sufficiently small and (4. 1 1) is obviously true. We need to investigate the case 

For e > we define an equivalence relation ~£ on by declaring p ~e q if and only if p and 
q belong to the same component of dC^. Similarly, we define an equivalence relation ~o on to 
declaring p ~o Q if and only p and q belong to the same connected component of dC^. Thus, for 
e > the number of connected components of dC^ is equal to the number of equivalence classes of 
~e. To prove the equality (4.1 1) it suffices to show that 



(4.12) 



Indeed, if (4.12) holds, then we deduce that the number of equivalence classes of ~£ is not larger 
than the number of equivalence classes of ~o- Since the number of connected components of Cq is 
independent of e if e is small and 

dist(C+,C+) ^0 as e\0 

we deduce that has at least as many components as Cq. 

Fix a component R of Cq and p,q & SR. We denote by [p, q]ji the arc of dR obtained by traveling 
counterclockwise from p to q. Along this arc there could be another crossing points Pq = p, ■ ■ ■ ,Pk = 
g G tg, an^ange in counterclockwise order. We set 

Pj ■■= -^eiPj)- 

Each of the arcs can is of one of the following two types 

I. A line segment contained in c- 

II. A sub-arc of 5C| g that intersects L^^c only at its endpoints. 




type II 



Figure 9. An elementary region S{I5, r) cut by a line L^ ^- Ttie intersection of ttiis 
region wittt ttie lower half-plane determined by L^ ^. has one component R whose 
boundary is decomposed in arcs of two types. 



If [pj_i,Pj]ji is of type I so that it is contained in c> then the points Pj^i and pj are also contained 
in L^^c and we denote by [Pj^i,Pj]R the oriented line segment going from Pj^i to pj. Cleai'ly pj^-^ ~ 



Pi 



Suppose now that [Pj^i,Pj]R is of type II. The points and divide the boundary dCe^e into 
two arcs, one of which approaches [Pj^i,Pj]R in the Hausdorff distance as e — )• 0. We denote this arc 
by [pj_i,Pj]ji. Lemma 4.8 implies that the arc [Pj_i,Pj]R intersects ^ only at its end points if e is 



PIXELATIONS OF PLANAR SEMIALGEBRAIC SETS 27 

sufficiently small. For such e's the arc [Pj_i,Pj]R lies on the same side of L^^c as \Pj-i,Pj]R so that 
Pjii Pj- By transitivity we now deduce that 

■^eip)=Po--ePl = ^s{q)- 

This proves (4.12) and thus proves (4.9) in the case when the elementary set Ce,o is nondegenerate. 

Case 2. The elementary set Ce,o is degenerate, i.e., /5e o = '^e,o- We denote by the set consisting 
of the endpoints of the graph of /3e o and the intersection of this graph with L^ c- Similarly, denote by 
(resp. J^) the set consisting of the endpoints of the graph of fi^^e (resp. re,e) and the intersection 
of this graph with the line c- As in Case 1 we can invoke Lemma 4.8 to obtain bijections 

We continue to use the notations Cq and C^^ introduced in the proof of Case 1. In this case is a 
finite union of subarcs of the graph of /5e o- Let these arcs be ^i, . . . , A}^. Each of these arcs carry a 
natural orientation. Denote by p^ the initial point of Aj and by the final point of Aj. 
We set 

p^j{e) := ^^(p,), Pl(^) := %iP,)- 

We define (e) and in a similar fashion. Consider the simple closed curve Aj which is the 

union of the following four aixs; see Figure 10 

• The line segment from Pj{e) to Pj (e)- 

• The arc of /3e,e from (e) to Qji^)- 

• the line segment from Qjis) to (e). 

• The arc of Te,e from (e) to pj(e). 




4- 

Figure 10. The arc Aj and the simple closed curve A5j, p'jis) — > Pji^) — ^ 9j (^) ~^ 
qr(e)^p-j{e). 

Lemma 4.8 implies that for e sufficiently small the closed curve Aj is contained entirely in the 
half-plane > c} so the bounded region it sunounds is contained in this half-plane as well. The 
region consists precisely of the regions surrounded by the closed curves Aj, j = 1, . . . , A; so that 
x(C'o^) = x{C^) = k for all e sufficiently small. This proves (4.9) in Case 2. 
Case 3» Cf> ^0 is a mixed elementary set. It has a minimal good partition 

a = to < ti < ■ ■ ■ < tn = h, 71 > 2, 

where for each j = 1, . . . , n the intersection of Ce,o with the strip [tj^i,tj] x M is either degenerate 
or nondegenerate. The intersection of the graphs of /? and r with each of the vertical lines x = tj, 
j = 0, . . . , n, is a singular point of Ce,o; see Figure 11. Since the c is not a critical value of the 
restriction of ^ to Ce,o, we denote that the line c does not contain any of these singular points. 
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Figure 1 1 . The minimal good partition of a mixed elementary sets. 

Set H^^ := > c}. For j = 1, . . . ,n and e > we set 

Ri := {[tj-i,tj] X M) n Clo, Ri,e := {[tj~i,tj] x M) n Ce,,. 
For k = 1, . . . ,n — 1 and e > we set 

Vk := {x = tfc} n Clo, Vk,e ■■= {x = tk} n C7e,e. 
To prove (4.9) it suffices to show that 

x{RjnH+J = x{Rj,enHlJ, Ve«l, j = l,...,n, (4.13a) 

x{VknH+J = xiVk,enH+J, Ve<l, A; = l,...,n-1, (4.13b) 

The equalities (4.13a) follow from the Cases 1 and 2 investigated above. The equalities (4.13b) are 
consequences of the following simple facts. 

• For any k = 1, . . . ,n — 1, the set Vk consists of a single point that does not line on the line 

• For any k = 1, . . . ,7i — 1, and e > the set Vk^s consists of a single vertical line segment. 

• For any A; = 1, . . . , n — 1, the set Vk^s converges in the Hausdorff metric to the set Vk. In 
particular, for e ^ 1 we have 

Vk C ^ Vk,e C 
This completes the proof of Theorem 4.4. 

□ 

Proof of Lemma 4.8. The inclusion (al) follows from the fact that the distance between the graph 
of Pe,e and the graph of /3e o approaches zero as e — )• 0. To prove (a2) let us denote by {xo,yQ) the 
coordinates of p. From the choice of r we deduce that for e sufficiently small the interval 

Je := [xq - r - a{e)e, xq + r + cr(e)e] 

is contained entirely in an interval of the form (cj_i, Cj) for some j = 1, . . . ,7i (where cj were defined 
much earlier in the proof to be the x-coordinates such that either /S^ or fail to be real analytic) so that 
/3e,o is on Jg. We set 

Ki = sup \/3'efi{x)\, K2 = sup \li%{x)\. 
We denote by the slope of ^ and by mo the slope of the tangent to the graph of /3e,o at p, 

Because ^ intersects the graph of /3e,o transversally at p we deduce ttiq 7^ m^. We deduce that for 
every x G we have 

IP'eflix) - mol < K2\x - xo\. (4.14) 
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The function /3e,e is piecewise linear. Consider the portion of this graph 
with the property that Pg, . . . ,p\^^^^_-^ are successive vertices on the graph of /3e ^ such that 

(pf , . . . C S^(p), pl.Pi^e) ^r{p)- 

Denote by (x^, y^) the coordinates of p^, j = 0, . . . , l{e), and set = /3e,o(3^j)- Observe that 

< (Ki+4)e, Vj = 0,...,^(e). 

We deduce that 

:= — = — 1- O ' 



where the constant implied by the O-symbol is independent of e. The mean value theorem implies that 
the difference quotient in the right-hand side of the above equality is equal to the derivative of /3e,o ate 
a point rf^ G (x^^^, xp. Using (4.14) we deduce that 

- mo I =0 + l^i-i ~ + - . (4.15) 

Since cr(e) — )• 00 we deduce that given 

7 < min r, ^|m^ — mo| 

there exist constants 6 = ^(7) > and £(7) > such that, for any e < £(7) the segments of the graph 
of /?£,£ situated in the strip |x — xq\ < 7 have slopes located in the range (mo — 7, mo + 7). In 
particular, none of these slopes can be equal to m^, and they are all situated on the same side of m^ as 
mo. 

If we fix 7 as above we can find ei (7) > such that, for e < ei (7) all the intersection points of L^^c 
with the graph of /3e ^ located in Sr(p) are in fact located in the nan^ow strip |j; — rro] < 7. The above 
discussion then shows that all these intersections must be transversal and they all have the same sign, 
sign(p, /3e,o)- Denote by Ne{'y) the number of such intersections. Set 

= ( Xo ± 7, /3e,o{xo ± 7) ) , = ( Xo ± 7, /3e,e{xo ± 7) ) . 

Consider now the closed curve obtained as follows. 

• Travel from the point P~ to P+ along the graph of /3e,e . 

• Next, travel on the vertical segment connecting P+ to Pq. 

• Travel along the graph of (3^,0 from P^ to P^. 

• Finally, travel along the vertical segment connecting Pq to P~. 

The above discussion shows that the intersection number between the line ^ and the curve 
is lb ( A'e (7) — 1 ) . On the other hand, since this curve is homologous to a circle, we have that the 
intersection number c and is 0. Therefore we conclude that N^i'j) = 1, which completes the 
lemma in the case of the function Pe,e- 

The above proof can be repeteated replacing /3e,£ with Te,e for the upper boundary case. □ 
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Appendix A. Semialgebraic geometry 
A set X C M" is called semialgebraic if it can be written as a finite union 

X = XiU-- -UXn, 

where each of the sets Xi is described by a finite system of polynomial inequalities. 

A map F : Xq — Xi between two semialgebraic sets Xi G i = 0, 1, is called semialgebraic if 
its graph Fj? is a semialgebraic subset of R"""*""!. 

Here is a Ust of basic properties of semialgebraic sets and functions. For proofs and more details we 
refer to [2, 3, 13]. 

• The union, the intersection and the Cartesian product of two semialgebraic sets are semialge- 
braic. 

• If X,Y are semialgebraic subsets of then so is their difference. 

• A subset of M is semialgebraic if and only if it is a finite union of open interval and points. 

• (Tarski-Seidenberg) The image and preimage of a semialgebraic set via a semialgebraic map 
are semialgebraic sets. 

• If / is an interval of the real axis and /:/—)• M is semialgebraic, then there exists a finite 
subset F C I such that the restriction of F to any component of / \ F is monotone and real 
analytic. 

• (Curve selection) If X is a semialgebraic subset of M" and xq € cl{X) \ X, then there exists 
a continuous semialgebraic map 7 : (0, 1) — >■ X such that 

lim 7(t) = xn. 
t\o 

• (Lojasewicz' inequality) Suppose that X is a compact semialgebraic set and /, : X — >■ M are 
continuous semialgebraic functions such that 

{f = 0}c{g = 0}. 

Then there exists a positive integer and a positive real number C such that 

|5(^)r <C|/(x)|, VxGX 

• Suppose that X is a compact semialgebraic set and / : X — > E is a continuous semialgebraic 
function. Then the function M — )• Z that associates to each t £ M the Euler chai-acteristic of the 
level set {/ = i} is a semialgebraic function. 

• A semialgebraic set is connected if and only if it is path connected. 

• A semialgebraic set has finitely many connected components and each of them is also a semi- 
algebraic set. 

Proof of Proposition 2.8 We prove only the statement about the total curvature. The statement about 
the perimeter follows the same pattern and has fewer complications. First some terminology. 
A continuous function / : [a, 6] — M is said to be piecewise if there exists a finite set 

S = {a = So < si < ■ ■ ■ < Si = , 

such that for any j = 1, . . . , i, and any k = 1,2 the restriction of / to the open interval Sj) is a 

function and the hmits 

lim f^^\x), lim f^^\x) 

exist and are finite. 

We say that the arc C is convenient if there exists a piecewise C^-function / : [a, 6] — > M such that 
either C is a graph of the function 

y = fix) 
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or C is the graph of the function x = f{y). When C is convenient, Proposition 2.8 is a special case of 
[12, Prop. 3.6]. 

To deal with the general case let us observe that since C is semialgebraic there exists an ordered 
samphng of C 

Q = {Qo,---,Qn} 

with the following properties. 

(a) The arc C starts at Qq and ends atQj^. 

(b) The arc C is smooth at each of the points Qi, . . . , Qn-i- 

(c) For any j = 1, . . . , A^, the portion of C between Qj-i and Qj is convenient. We denote by Cj 
this portion. 

Denote by the ordered sampling of Cj determined by the points in contained in Cj. We 
denote by Ki the total curvature of the PL-curve C{y^). Since each of the curves Cj is convenient 
we have 

\miKi = K{Cj), Vi = l,...,iV, 

e\0 

SO that 

N 



\unJ2Ki=K{C). 



On the other hand, since C is smooth at the points Qi,...,(5Ar_iwe deduce that 



Appendix B. The approximation algorithm 



□ 



In this section we give a more formal description of the approximation algorithm. We convert Pe{S) 
into an m X m matrix A (where m depends on e) of I's and O's, where = 1 if and only if the 

pixel of center Cij{e) touches S. 

Given this matrix A we will generate a PL set Se which approximates the original set 5. We will 
assume that e is fixed thi^oughout the description of the algorithm. 

The algorithm depends on two parameters, both positive integers: the spread a and the noise width 
V. These should be chosen using functions a{e) and v{e) which have the properties: 

lim ea{e) = 0, lim e{a{e))'^ = oo 

lim eu(e) = 0, lim iy(e)e^~''° = oo 

where kq is a constant dependent on S introduced in Theorem 3.4. However, for purpose of performing 
the algorithm it is only relevant that we have chosen these two integers a and u. 

The algorithm uses several smaller subroutines. The first subroutine Stack obtains information about 
the various columns of A which will be used to determine both the noise intervals as well as to select 
the vertices of Se- The input of Stack is a list 

C = Ci, . . . Cm, Ci = 0, 1, 

where C is one of the columns of A. The output of stack is a list of nonnegative integers 

n(C); 6i < ti < 62 < ^2 < • • • < K{C) < tn(C), 
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where n{C) is the number of stacks in the column encoded by C, and the location of the bottom and 
top pixel in the j-th stack is determined by the integers bj,tj. More formally 

Cfc = 1 31 < j < n(C) : bj < k < tj. 

If C = Ci, the i-th column of A, i.e., 

Ci — ) ■ ■ ■ 5 CLi,mi 

then we will denote the output Stack (Cj) by 

A number l<i<m — lis called a jump point if 

rij / rij+i. 

The next subroutine is called jump. Its input is an integer k e [1, m) and the output is an integer 
jk = \urr\p{k) where is the next jump point, i.e., 

|« G [k, m) n Z; i is a jump point } = 0, 

then we set 

jump(/c) := m + 1. 

Otherwise 

jump(A;) = minji G [A;, m) n Z; i is a jump point }. 

Using these subroutines we can construct the noise regions of the approximation. These are simply 
the columns which are within 2u columns of a jump point. Specifically we create a certain number of 
intervals: 

[ii,ri], [ia,ra] C [l,m] 
where the integers £k, fk are determined inductively as follows. 

h = max(jump(l) - 2u,l), ] ri = min( m, jump(l) + 2iy). 
Suppose that ii,ri, . . . , £j,rj are determined. If jump(rj) > m we stop. Otherwise we set 

£j+i = max(jump(rj) -2u,l), rjj^i = min( m, jump(rj) + 2u). 
The intervals [ii,ri], . . . , [la, may not be disjoint, but their union is a disjoint union of intervals 

[ai,6i], . . . , [aj,6j], h < Oj+i. 
The intervals [aj, 1 < j < J are the noise intervals. The intervals 

[1 , ai] , [6i , a2] , . . . , [bj-i, aj] , [bj, m] 

are the regular intervals. 

Now that we have determined the noise and regular intervals, we can create the approximation Sg. 
We do this with separate procedures on the noise or regular intervals. In either case the approximation 
will be formed by (possibly degenerate) trapezoids whose bases are vertical.We call any set which is a 
union of finitely many such trapezoids a polytrapezoid. The approximations on the regular and noise 
intervals will both be poly trapezoids, and Sg itself will also be a polytrapezoid. 

First some notation. Given a collection of points 

Bo,T(j,. . . , Bn, T/v G 

such that 

x{Bi) = x{T,), y{Bi)<y{T,), Vi = 0,...,iV, 
x{Bj^i) < x{Bj), VI < i < N, 
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we denote by polygon(i3o, To, . . . ,Bn,Tn) the region sun-ounded by the simple closed PL-curve 
obtained as the union of line segments 

[Bo, Bi],. . . , [Bn~i,Bn], 

[Bn,Tn], • • • , [Ti,To], [To,Bo]. 
Note that each of the quadrilaterals Bi^iBiTiTi^i is a (possibly degenerate) trapezoid with vertical 
bases. 

Consider first the regular intervals. Given a regular interval / := [p, q] we observe that the number 
of stacks rii is independent of i G [p, q]. We denote this shared number by n = n{I). 
We construct inductively a sequence of numbers < • • • < as follows: 

• We set iQ = p. 

• If (/ — p < 2(7 we set = 1 and ii = q. 

• If iQ, . . . ,ik are already constructed, then, if q — ik < 2a we set N = k + 1 and ik+i = q, else 
ik+i = ik + (y- 

Note that if q — p> a, then N > \, = p, ij^ = q and 

A^ = l if q - p < a. 

We have 

stack(Cjj.) = n, . . . , &jj.,n) Uk,n- 

For j = 1, . . . , n, and A; = 0, . . . , we denote by Bi^j the center of the e-pixel corresponding to the 
element entry bi^,j in the column Cj^, . Similarly we denote by T^j the center of the pixel corresponding 
to the entry ti^j of the column Ci^. For 1 < j < «,(/), we set 

■Pj{I) := polygon (Bo,j, 7o,j, • • • , Bnj^Tnj). 

Define 

ni) = U ^jW' ■■= U ^(^)- 

j=l I regular interval 

Suppose now that / = [p, g] is a noise interval. We modify the column 

Cp — 0,p,l, • • • ; Ojp,m 

to a column 



Cp — CLp ]^ , . . . , ffl. 



/ 

• • • ) '^p^mi 



by setting 

0'p,k '■= ' 

We apply the subroutine stack to the new column C^and the output is 

Stack(Cp) = n(/), bi <ti < ■ ■ ■ < bn <tn- 

For j = 1, . . . , n{I) we set 

Bo J := A[p,bj], Toj := A\p,tj], 

Bi,j ■■= A[q,bj], Toj := A[q,tj], 
(recall that ^[«, j] is defined as the center of the pixel associated to (aij). Next, for j = 1, . . . , nil) 
we define the rectangle 

:= polygon(Po,i,7bj,5i,i,Tij), 
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U ^w- 



j=i I 

The output of the algorithm is the polytrapezoid 



/ noise interval 




noise • 
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